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o ■ 

' In a previous paper [Phys. Rev. A 72, 033415 (2005)], it was shown that sub-Doppler coohng 

^SJ , occurs in a standing-wave Raman scheme (SWRS) that can lead to reduced period optical lattices. 

These calculations are extended to allow for non-zero detuning of the Raman transitions. New phys- 
ical phenomena are encountered, including cooling to non-zero velocities, combinations of Sisyphus 
and "corkscrew" polarization cooling, and somewhat unusual origins of the friction force. The cal- 
. culations are carried out in a semi-classical approximation and a dressed state picture is introduced 

' to aid in the interpretation of the results. 

^ PACS numbers: 32.80.Pj,32.80.Lg,32.80.-t 
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^ ■ I. INTRODUCTION 

o : 

' In a previous paper I*] (hereafter referred to as I), we have shown that sub-Doppler coohng occurs in a standing-wave 
^ Raman scheme (SWRS). The SWRS is particularly interesting since it is an atom- field geometry that leads to optical 
O ■ lattices having reduced periodicity. Reduced period optical lattices have potential applications in nanolithography 
c/) [ and as efficient scatterers of soft x-rays. Moreover, they could be used to increase the density of Bose condensates 
J>~». in a Mott insulator phase when there is exactly one atom per lattice site. With the decreased separation between 

■ lattice sites, electric and/or magnetic dipole interactions are increased, allowing one to more easily carry out the 
entanglement needed in quantum information applications 1^. 

In this paper, the calculations of I, which were restricted to two-photon resonance of the Raman fields, are extended 
T-H to allow for non-zero Raman detunings. There are several reasons to consider non-zero detunings. From a fundamental 
^ ' physics viewpoint, many new effects arise. For example, one finds that, for non-zero detuning, it is possible to cool 
I atoms to non-zero velocities, but only if both pairs of Raman fields in the SWRS are present, despite the fact that 
the major contribution to the friction force comes from atoms that are resonant with a single pair of fields. This 
is a rather surprising result since it is the only case we know of where non-resonant atoms that act as a catalyst 
for the cooling. Moreover, comparable cooling to zero velocity and non-zero velocities can occur simultaneously, but 
\^ the cooling mechanisms differ. We also find effects which are strangely reminiscent of normal Doppler cooling, even 
, though conventional Doppler cooling is totally neglected in this work. A dressed atom picture is introduced to simplify 
' the calculations in certain limits; however, in contrast to conventional theories of laser cooling, nonadiabatic coupling 
O ' between the dressed states limits the usefulness of this approach. The non-adiabatic transitions result from the unique 
'U^ potentials that are encountered in the SWRS. To our knowledge, there are no analogous calculations of laser cooling 
in the literature. 

1 From a practical point of view, there is also a need for calculations involving non-zero detunings. For example, in 
the quantum computing scheme proposed in 0] , the Raman frequency differs at different sites owing to the presence 
^ : of an inhomogeneous magnetic field, making it impossible to be in two-photon resonance throughout the sample. As 

a result, one has to assess the modifications in cooling (and eventually trapping) resulting from non-zero detunings. 
^ , The basic geometry is indicated schematically in Fig. 1. Transitions between states |1) and |2) in the Raman 

■ scheme occur through the common state |3) using two field modes. Consider first the effect of fields Ei and £2- Field 
El, having frequency f2i and wave vector ki = k =kz drives the 1 — 3 transition while field E2, having frequency 
fl2 ~ ^1 — LO21 and wave vector k2 ~ — k drives the 2 — 3 transition, where iOji is the frequency separation of levels 
j and i (it is assumed that Q.2/C ~ i^i/c, or, equivalently, that W21/W31 <C 1). Owing to polarization selection rules 
or to the fact that 0)21 is greater than the detuning A = fii — wjj, one can neglect any effects related to field Ei 
driving the 2 — 3 transition or field E2 driving the 1 — 3 transition 3] . If, in addition, the atom-field detunings on the 
electronic state transitions are sufficiently large to enable one to adiabatically eliminate state |3), one arrives at an 
effective two- level system in which states |1) and |2) are coupled by a two-photon "Raman field" having propagation 
vector 2k and two-photon detuning 5 = — Q.2 — ^21- 

Imagine that we start in state |1). If the initial state |1) amplitude is spatially homogeneous, then, after a two- 
quantum transition, the final state ( state |2)) amplitude varies as e'^*'^^. Such a state amplitude amplitude does 
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FIG. 1: Schematic representaion of the energy level diagram and atom - field interaction for the standing wave Raman scheme 
(SWRS). 

not correspond to a state |2) population grating, since the final state density is spatially homogeneous. To obtain a 
density grating one can add another pair of counter-propagating fields as shown in Fig. 1. These fields E3 and E4 
differ in frequency from the initial pair, but the combined two-photon frequencies are equal, 

S ^ Cli - D,2 - UJ21 ^ ^3 - ~ W21. (1) 

The propagation vectors are chosen such that ks ~ — k4 = — k. The frequencies of fields Ei and E3 are taken to be 
nearly equal, as are the frequencies of fields E2 and -E4, but it is assumed that the frequency differences are sufficient 
to ensure that fields Ei and E3 (or E2 and E4) do not interfere in driving single photon transitions, nor do fields Ei 
and E4 (or E2 and E3) drive Raman transitions between levels 1 and 2 4]. On the other hand, the combined pairs 
of counter-propagating fields {Ei and E2) and {E^ and E4) do interfere in driving the 1 — 2 Raman transition and 
act as a ^^standing wave^^ Raman field which, to lowest order in the field strengths, leads to a modulation of the final 
state population given by cos(4A:z). In this manner, a grating having period A/4 is created. 

The friction force and diffusion coefficients are calculated using a semiclassical approach. For (5 7^ 0, they differ 
qualitatively from the corresponding quantities obtained in standard Sisyphus cooling. The physical origin of the 
friction force was discussed in I. The calculation can also be carried out using a quantum Monte-Carlo approach, but 
the results of such a calculation are deferred to a future planned publication. 

II. SEMI-CLASSICAL EQUATIONS 

As in I, we consider the somewhat unphysical level scheme in which states |1) and |2) in Fig. 1 have angular 
momentum J = 0, while state |3) has angular momentum J = 1. The field intensities are adjusted such that the Rabi 
frequencies x (assumed real) associated with all the atom-field transitions are equal (Rabi frequencies are defined 
by —/ii?/2fi,, where ^ is a component of the dipole moment matrix element between ground and excited states), and 
the partial decay rate of level 3 to each of levels 1 and 2 is taken equal to r/2 (equal branching ratios for the two 
transitions). The fields all are assumed to have the same linear polarization; there is no polarization gradient. The 
results would be unchanged if the fields were all a+ polarized. 

It is assumed that the electronic state detunings are sufficiently large to satisfy 

rii - UJ31 « f^a - UJ31 « - (^32 « - t^32 = A > r, X- (2) 

In this limit and in the rotating- wave approximation, it is possible to adiabatically eliminate state |3) and to obtain 
equations of motion for ground state density matrix elements. With the same approximations used in I, one obtains 
steady-state equations, including effects related to atomic momentum diffusion resulting from stimulated emission 
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and absorption, and spontaneous emission. In a field interaction representation , the appropriate equations are 
d{p22 - pii) 



a- 



dx 



dpi2 



dx 

P21 = P21 
or, in terms of real variables, 
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dx 

where the total population S evolves as 

dS 7 
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du 



= -h^k'^V —- ~ AV'ahksmx- h^k'^V cos x— 



dt 



and 



dp^ dp 



u = P12+ P21 , 
V = i {p2i - P12) , 
w = P22- Pii , 
S = P11+ P22 , 



with 



d = A 

a = fcv/r' 
a = A/r, 



r' 



A2 + (r/2)' --x'r/A' 



(3a) 

(3b) 
(3c) 

(4a) 
(5) 



(6a) 
(6b) 
(6c) 
(6d) 



(7a) 

(7b) 

(7c) 
(7d) 

(7e) 



Each of the functions u, v, w, S are now functions of the z-component of momentum p = My (Af is the atom's mass 
and V is the z-component of atomic velocity) as well as x, but it is assumed in this semiclassical approach that S is 
position independent. The parameter cr = A/F is assumed to be large compared with unity. 

It will also prove useful to define dimensionless frequencies normalized to Ur, momenta normalized to hk, and 
energies normalized to fiwr, where ujr is the recoil frequency 



= hk'^/2M, 

such that 6 = 6/oJr, f = T/uJr, f' = T'/uJr, Xeff = Xeff/<^r {Xeff = X^^/ 

two-photon Rabi frequency}, p — p/hk, etc. In terms of these quantities, 

X^^/uJr 



A2 + (r/2)^ 



(8) 



= rV) is an effective 



Xeff 



I, 



A2 + iT/2Y 
f' = I/a, 
A ^ la 121 , 
a = p/f = 2apll . 

Note that / is the effective coupling strength normalized to the recoil frequency. 
Equation (|4a|l can be written in matrix form as 



dB(a;) 

dx 



[Ai + 2(7 cos a; A2] B(a;) - A(a;), 



(9b) 
(9c) 
(9d) 



(10) 
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where 



B(a;) 





Aix) 



, Ai 



cos X S{p) + 2ahk sin 



' dp 




(11) 



(12) 



It should be noted that Eq. (|10|) differs qualitatively from the corresponding equation encountered in liigh intensity 
laser theory. Owing to the fact that decay of u, v, w is linked to spontaneous emission, the decay parameters depend 
on field intensity. When all frequencies are normalized to the optical pumping rate F', the effective coupling strength 
a is actually independent of field strength; moreover, since it is assumed that ct > 1, one is always in a "high intensity" 
limit. In contrast to the equations describing conventional Sisyphus cooling or high intensity laser theory, there is a 
source term for u, but no source term for the population difference w. 

The formal solution of Eq. H10|l satisfying boundary conditions resulting in a periodic solution is 



B(x) = -- / cia;'e-^i(^-="')/" [2acosa;'A2B(x') + A{x')] , 
a J-oo 

which, in terms of components, can be written as 

/•OO 

u{x) = dre^^ {2(TC0s(ar — x) sin [dr] w{x — ar) + s{x — ar) cos [dr]} 

Jo 

/>OC 

v{x) = dre^'^ {2crcos(ar — x) cos [dr] w{x — ar) — s{x — ar) sin [dr]} 

Jo 



where 



w{x) — 2<T / dre cos{aT — x)v{x — ar), 
Jo 



dS 

s(x) = cos X S(p) + 2ahk sin x— — 

op 



(13) 

(14a) 
(14b) 
(14c) 

(15) 



Substituting v{x) into the equation for w{x) we obtain 

w{x) = — 4ct^ / dre^^ cos(Q;r — x) dr' e^'^ cos(q;(t + t') — x) cos [dr'] w[x — a[T + r')] 
Jo Jo 

f^OO /"OO 

+ 2(T / dTe~'^ cos(aT — x) dr' e~'^ s[x — a(r + r')] sin [dr'] . 
Jo Jo 

Once the solution for w{x) is obtained, it is substituted into Eq. H14a(l to determine u(x) and the solution for u(x) 
substituted into Eq. (0 for dS/ dt. The resultant equation is averaged over a wavelength resulting in 



(16) 



= h k r - 4r ahk—r]i - -h k T -^m ; 
ot 5 op'^ op 5 ap^ 



where 



1 r^- 

771 = u{x) sinx = — / dxu{x) sin 
27r Jo 



r]2 = u{x) cos x = — / dxu{x) cos x . 
27r Jo 



(17) 

(18a) 
(18b) 



and the bar indicates a spatial average {S — S, by assumption). Equation l|17|) is then compared with the FokkerPlanck 
equation 



dS_ _ d_ 

dt dp 



8S 8 



(19) 



to extract the spatially averaged friction stimulated diffusion D^t, and spontaneous diffusion Dsp coefficients. 
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III. SOLUTIONS 



A. Backward recursion method 

Equation (|16|) can be solved using Fourier series and a backwards recursion scheme 0, Q, as outlined in the 
Appendix. In this manner one obtains 

F = ~2hkaT'if , (20a) 

Dst = ih^k^a^r'ist , (20b) 

Dsp = h'k^r' Q + A^^p^ , (20c) 

where £^sp, £,st are given as Eq. (A17) in the Appendix. An analytic solution for the ^s can be found only for 
d = [£,f{d = 0) = a(l + Q!^)"^ , £,sp{0) — Cst(O) = (1 + o;^)""^] ; otherwise, these quantities must be obtained via 
the recursive solutions. The effective field strength parameter in this problem is a and one might expect that a /a 
recursions are needed to solve Eqs. IjAGfl accurately ]EI3- Actually, the number of recursions required depends in a 
somewhat complicated manner on the values of several parameters. Each recursion introduces resonances at specific 
values of a which can be interpreted as Stark-shifted, velocity tuned resonances Q. For example, with d = 0, the 
lowest order recursive solution has a very strong (proportional to a^), narrow resonance at — 3/5, but the second 
order approximation removes the divergence, while introducing yet a second resonance. Some of these velocity tuned 
resonances are seen in some of the graphs presented below. For d > 1, an upper bound for the number of recursions 
required to map out all the resonances is of order (cr/d) (cr/a) ; for d<Clord/(T3>l only a few terms are needed. 
Even if only a single recursion is needed, the general expressions for ^sp, S,st are still fairly complicated. 

For d <C 1 one finds corrections of order d^ to the analytical results, but owing to their complexity, these expressions 
are not given here. For d/a 3> 1, one finds that, near a = 

2q! /ctn2 

^^-Tw(d) ' 

£,st 



1 / <j\ 2 

3 + 8a2) 
l + 4a2 



(3 + 8^. ay 



and near a = ±d 



2(l + /|+a2 



- .7.^^,, ^ (22a) 



^st - , i , ,^ , (22b) 



2(1 + /I 
1 + ^ 
2(l + /|+a2 



7777^:^7^^^ (22c) 



where f± — (a =F d). 

In the limit d ^ cr > 1, the friction force as a function of a contains three dispersion-like structures centered at 
a = 0, ±d. This implies that atoms can be cooled to these values of a. The amplitude of the component centered 
at a = is of order cr^/d^ while its width is of order unity. On the other hand, the amplitude of the components 
centered at a = ±d are of order 1/a while their width are of order a. It is shown below that the central peak in the 
momentum distribution is negligible compared with the two side peaks in the limit d 3> cr > 1; that is, in this limit 
cooling occurs more efficiently to velocities a = ±d for which the atoms are Doppler shifted into resonance with the 
two-photon transitions connecting levels 1 and 2. The width of the a = ±d components is similar to that found in 
sub-Doppler cooling in magnetic fields (MILC) 0; in both MILC and the SWRS, the qualitative dependence for the 
friction coefficient as a function of a is similar to that found in sub-Doppler cooling using " corkscrew" polarization 
[l0|. As such, the curve is "power broadened," since the effective field strength in the problem is cr. 

It is tempting to consider the contribution to the friction force near a = d as arising from the single pair of fields 
(El and E2), since these fields are nearly resonant with the 1-2 transition in a reference frame moving at 2fcv = S. 
Tempting as it may be, this interpretation is wrong, since we have already shown in I that, for a single pair of Raman 
fields, the friction force vanishes identically, regardless of detuning. Thus, it is necessary that the second pair of fields 
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FIG. 2: Averaged force and diffusion coefficients in dimensionless units as a function of scaled momentum p = p/hk for / = 10, 
a = 10, and S = 20. 



be present, even if they are far off resonance with atoms satisfying 2fcv = S. The main effect of the second pair of 
fields is to cancel the contribution to the force from the population difference between levels 1 and 2 (see Appendix 
A in I), leaving the contribution from the 1-2 coherence only (u = pi2 + P2i)- Near 2A:v = 5 the major contribution 
to u does come from atoms that are nearly resonant with the 1-2 transition in a reference frame moving at 2/cv = 5, 
but the scattering of the second pair of fields (-E3 and E4) from the population difference created by the resonant pair 
of fields modifies the net force on the atoms. In some sense, one can view the second pair of fields as enabling the 
cooling at 2/cv = 5. Note that the magnitude of the damping coefficient is down by from that at d = 0; it is of the 
same order as that found in sub-Doppler cooling using " corkscrew" polarization jl^l . 

For arbitrary values of a and 5, with a of order 10, the recursive solution converges very rapidly for most values 
of a and numerical solutions can be obtained quickly and easily. Two examples are shown in Figs. 2 and 3, where 
the averaged friction force F in units of hkoj^ and the averaged diffusion cocfhcients Dgt and Dgp in units h^k'^uj^ are 
plotted as a function of the scaled momentum p = p/hk. In terms of the ^s defined by Eqs (|20l) . these quantities can 
be written as 

F/hkujj. 

Dst/h^k^UJr 
Dsplh^k^UJr 

In Fig. 2, / = 10, a = 10, and 5 = S/ujr = 20. One sees in these curves a number of velocity tuned resonances under 
a "power-broadened" envelope jg. In Fig. 3, / = 10, ct = 10, and 5 — 220, implying that d = 110 and a/d = 1/11. In 
this limit Eqs. 1)21(1 . 1)22(1 are valid and we see three contributions to the averaged force and diffusion coefficients. The 
values of the force and diffusion coefficients near the Doppler tuned resonances aX p — ±5/4 are typical of corkscrew 
polarization cooling "T^, and the ratio of the force to diffusion coefficient is of order 1/hk. On the other hand, this 
ratio is of order l/hka near p = 0, a result that is typical of Sisyphus cooling in a lin_Llin geometry; however, both 
the friction and diffusion coefficients are smaller than those in conventional Sisyphus cooling by a factor (cr/d)^ when 
a/d ^ 1. As a consequence, the cooling is dominated by the contributions near p = ±5/4 when cr/d ^ 1. 
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FIG. 3: Same as Fig. 2, with / = 10, cr = 10, and S = 220. 



B. Iterative Solution 



Since the effective field strength is always greater than unity, perturbative solutions of Eqs (|l()|l are not of much use. 
However, one can get a very rough qualitative estimate of the dependence on detuning of the friction and diffusion 
coefficients near a = by considering an iterative solution of Eqs. (|10|l in powers of a. This will work only in the 
limit that a < 1, so it cannot correctly reproduce the contributions to the friction and diffusion coefficients at a = ±d 
when d 3> cr > 1. The iterative solution is useful mainly when d^ < cr^, since, in this limit, the dominant contribution 
to the momentum distribution comes from the region near a = 0. 

The iterative solution is straightforward, but algebraically ugly. To order a, one obtains from Eqs. H10|l 



where A(x) = [Ai + 2(tcosx A2]. When the u component of B(a;) is extracted from this solution and the result is 
substituted into Eqs. (|18|l . all the integrals can be carried out analytically and one finds 




(23) 
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where 

n = 9.rrfik 

dp 



dS 

2ahk— , (24a) 



- 1 - . , (24b) 



7d{ld+la) 



7<T (7d + 7<t) 



2d 



2 



7;?72 



Id 



VT+d2 , 7^ = Vr+d2+4a2. (24e) 



By comparing Eqs. (|17|l . p9|) and neglecting the contribution from the second term in the equation for 772 (since it is 
of relative order hk/p), we extract the spatially averaged friction and diffusion coefficients 

F = -2hkaT' aS^i'^ , (25a) 
D^t = 4n2fc2^2p,^(0) ^ ^25b) 

Ds, - n^fc^r' f I + Ae^"^ • (25c) 
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These are all even functions of the detuning d. 

The spatially averaged form factors S,2^\ are equal to unity for d = 0, but vary as 

«ei°^«(3/2)d'^^3aVd^ (26) 

for d ^ d/(T ^ 1, in agreement with Eqs. 121() . In this limit, both F and Dind approach zero, but Dgp approaches a 
finite value since Rayleigh scattering of the fields is independent of 6 ioi S A. The friction force when d ^ d/cr ^ 1 
is given by 

F ~ -AT'hkaaa'^/d^ = -IQhk^v [r (A^F^)] {x^ / Af 

This equation is strangely reminiscent of the equation for Doppler cooling of two-level atoms by an off-resonant 
standing wave field where one finds 

FDc---^nk{kv)T{x^) /A\ 

taking into account the fact that twice the momentum is transferred in a two-photon process. For the expressions to 
agree, one must associate a "two-photon spontaneous scattering rate" Ftp — F (A(5/r^) with the Raman transitions. 

Of course, if d/cr ^ 1, the contributions to the friction and diffusion coefficients near a — ±d become dominant 
insofar as they affect the momentum distribution. In this limit one cannot use the iterative solution to estimate the 
equilibrium temperature since the contributions from higher velocity components play a significant role. 



C. Dressed State Solution 



The effective Hamiltonian for the SWRS, neglecting decay is 

jj ^ll [ ^ 4Xe//C0SX 

2 V 4Xe//cosa; -S 

By diagonalizing this Hamiltonian one obtains semiclassical dressed states whose energies, as a function of x are said 
to characterize the optical potentials associated with this problem. It turns out that the use of dressed states in the 
SWRS is of somewhat limited use, owing to nonadiabatic coupling between the potentials. Nevertheless, the dressed 
states do provide additional insight to the cooling dynamics. 
The eigenvalues of H are given by 




A± = ±hR/2, R = + 



(28) 
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FIG. 4: Dressed state potentials. 



along with eigenkets 



where 



1^) = c|l) +s|2) , (29a) 
\B)^-s\l)+c\2), (29b) 



c = cos9, s — sin9, (30a) 

Ro = 4xe/ / cos X, (30b) 

cos(26l) = 6/R, sin(26i) = Rq/R, O<0 < tt/2. (30c) 

The optical potentials are sketched in Fig. 4. As 6/2xeff tends towards zero, the potentials "touch" whenever 
cos a; — O. As is seen below, nonadiabatic transitions occur at such points 
Defining dressed state amplitudes via 

ao = Ta (31) 

with 

-^(::)^ ^^<z)- ^^-.i 



10 



and a dressed state density matrix po = a/^aj^, one can transform Eqs. © into the dressed basis as 



dwp 
dx 



dx 



as 

'dt 



~WD + sin(26l) 
sin(26')i;£i 



cos(a;)5 + 2ahksm {x) 

de 



dS_ 
dp 



2a 



2a— {pAB + pea) 



- (1 + iD) PAB - cos(26i) 



— cos{x)S + ahk sin (x) 



2a 



{PBB - Psa) , 
ox 



h^k^r'^ ahk sin. 



d [cos{20)ud - sin(26')w£,] 



dp 



h^k'^T' cos X 



(92 [cos(26')ui5 - sin(26l)wD] 



5 ' dp^ 
ud ^ PAB + Pba, vd ^ Hpba - Pab) , wd ^ pbb - Paa, 

PBA = P*AB ' 

D = R/2T'. 



(32a) 



(32b) 



(32c) 

(32d) 
(32e) 
(32f) 



For (7^1, the terms varying as a~^ can be dropped. If one also neglects the nonadiabatic coupling proportional to 
dO/dx, Eqs. I|32|l have the remarkable property that, even in the presence of dissipation, the equations for the dressed 
state coherences and populations are completely decoupled. Assuming for the moment that such an approximation is 
valid, one has the immediate solution 

1 dS 
- cos{x')S + ahk sin (x') — 



PAB — — {oi) ^ / da;' cos [26' (x')] cos(a;')S' + crfi.A: sin (x') ^ exp [— (1 + iD) (a; — x') /a] ; (33a) 

— oo 

dx' sm[29 {x')] cos{x')S + 2ahksm{x) — exp [— (x — a;') /a] . (33b) 

- oc L P . 

It then follows from Eqs. H32() . and H19() that the spatially averaged friction and diffusion coefficients are given by 



D 



F = 4fifccrr'sin {x) {cos [26* {x)\ {C{x) + C*(x)\ - sin [20 (x)] A(a;)}, 

-%h^k^a'^V'sm (x) {cos [261 (x)] [^(a;) + D*{x)\ - sin [20 (x)] B(x)}, 
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cos (x) {cos [26* (x)] [C{x) + C*(x)] - sin [26* (x)] A{x)} 



(34a) 
(34b) 

(34c) 



where 



A{x) = [a) ^ I dx' sin [20 (x')] cos(x') exp [— (x — x') /a] , 

/X 
dx' sin [20 (x')] sin(x') exp [- (x - x') /a] , 
-oo 

/X 
dx COS [20 (x')] cos(x') exp [- (1 + iD) (x - x) /a] , 
-oo 

/X 
dx' COS [20 (x')] sin(x') exp [- (1 + iD) (x - x') /a] , 
-CXj 

and the bar indicates a spatial average. In general, the integrals and spatial averages must be calculated numerically. 

In contrast to other dressed state theories, the dressed states here are of limited use since the nonadiabatic coupling 
is always significant. This is related to the fact that the decay constants are intimately related to the coupling strength, 
that the potentials periodically approach one another, and that the nonadiabatic coupling is maximal at these close 
separations [dO/dx ~ (c/d) sinx]. The dressed picture gives a reasonable approximation to the friction and diffusion 
coefficients when |a±d| ^ 1 and d > 3> 1. In this limit one can make a secular approximation and ignore the 
contribution from the C(x) and D{x) terms in Eqs. H34|) . The nonadiabatic terms neglected in Eq. (|32|l are of order 
aa/A^ in this limit. Thus, the approximation is valid for relatively large detunings and values of a less than or on 
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(b) 

(a) 



S(P) 




FIG. 5; The momentum distribution S{p) as a function of 5 = 5/ujr for / = 40, a — 10; 3-dimensionaI plot (a) and density plot 
(b). 

the order of unity. Indeed, in the limit cr/d ^ 1, the dressed picture results reproduce those of Eq. H21() provided a 
is not too large. On the other hand, they do not reproduce those of Eq. (|22|l near the Doppler shifted two-photon 

resonances; the dressed results vary as (l + /^) ^ rather than (l + + cr^) ^ . Both the secular approximation and 
the neglect of nonadiabatic coupling break down near these two-photon resonances. 

For the nonadiabatic coupling to be negligible compared with convective derivatives such as adwo/dx, it is necessary 
that do / dx ^ 1. It can be shown that in the regions of closest approach of the potentials that 89/ dx ~ Xeff/^ = cr/d. 
Thus, for the dressed picture to be valid, one is necessarily in the region where the approximate solutions Eas. (|21|l . 
(|^ are all that is needed. 



D. Density matrix solution 



As a final approximate approach one can adiabatically eliminate pi2 and p2i from Eqs.||2Jl. This procedure will 
allow one to obtain an analytical solution for all density matrix elements in terms of a sum over Bessel functions. 
Such an approach is valid for 6 ^ V and 5 l>d I ^ so it has a limited range of applicability. The detailed results 
are not presented here. 



IV. MOMENTUM AND ENERGY DISTRIBUTIONS 



In terms of the normalized momentum -p ~ p/hk, the steady state solution of the Fokker-Planck equation. 



hk dp 

subject to the boundary condition dS/dp\p=o = 0, is 



DtotdS f- 1 dPsp ^ ^ 



hk dp 



(35) 



S(p) = So exp <hk dp 



_i y hk dp' J 



D 



tot 



(36) 



where 



St) 



dp exp ■[ hk I dp' - 



hk dp' 







(37) 
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Taking into account definitions, Ea. H2U|l . we obtain 



S{p) = ^oexp ■ 



dp'^ o 

5 10^' 



(38) 



where the are define in Eqs.(A17) and we neglect the term {1/hk) dDgp/dp' . 

The momentum distribution S{p) is plotted in Fig. 5 as a function of 5 for cr = 10 and / = 40 and in Fig. 6 
as a function of / for cr = 10 and 5 — 220. The curves in Fig. 7 are cuts of Fig. 5 for 5 — 0, 170, 220, and 300. 

/ _ \l/2 

When d/(T = 5/21 <C 1 and / ^ cr > 1, there is a central component having width of order \/2T I 1 + P/8I^ ] , 

that is estimated using Eqs. (|24(l . For d/cr = 5/21 ^ 1, the momentum distribution breaks into three components 
centered at p = 0, ±(5/4 [2fcv = 0,± 5], with the central component negligibly small compared with the side peaks 
{relative strength of side to central peak scales roughly as (d/a)"'^^^'', estimated using Eqs. (|21() . (|22|l l. The width 
of the side peaks for d/cr = 6/21 ^ 1 also scale as \/7, although they are slightly broader than the central peak 
when d/cr — 5/21 ^ 1, reflecting the fact that the side peak cooling is of the corkscrew polarization nature, while 
the central component cooling for d/cr ^ 1 is of the Sisyphus nature. For intermediate values of d/a three peaks in 
the momentum distribution are seen clearly; for example, when 5 = 220, / — 40, a — 10, the amplitudes of the three 
peaks are equal. 

The mean equilibrium kinetic energy can be calculated according to 



E^E,g/Er^ / dpp^Sip) 



(39) 



where Er — Hujr is the recoil energy. This quantity must be calculated numerically, in general. However, for 

/ , \l/2 

d/cr = 5/21 < 1 and / > cr > 1, one can estimate that E = I [l + 5"^ /M] , using Eqs.lI23|). For d/cr = 5/21 > 1, 



the side peaks lead to an equilibrium energy that scales as since momentum components at both p — ±(5/4 

are present; however, the energy width associated with each side peak scales as /. In Fig. 8, we plot E — Eeq/Er as a 
function of / for cr = 10 and several values of 5. 
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FIG. 7: Cuts in the momentum distribution of Fig. 5 for S = O(dash-dotted line), 170(dotted line), 220(daslied line), 300(solid 
line) . 



V. SUMMARY 

We have extended the calculations of I to allow for non-zero detuning in a standing-wave Raman scheme (SWRS) 
that results in reduced period optical lattices. The results differ from that of conventional Sisyphus cooling. Optimal 
cooling occurs for exact two-photon resonance, but many new and interesting features appear for non-zero detuning. 
A dressed atom picture was introduced, but had limited usefulness, owing to the presence of nonadiabatic transitions. 
In a future planned publication, we will look at Monte Carlo solutions to this problem and examine the dynamics of 
the velocity distribution. Specifically we will attempt to determine how the atomic momentum jumps between the 
momentum peaks shown in Fig. 5. Furthermore we will see if it is possible to localize atoms in the potential wells 
shown in Fig. 4. The ability to do so would imply separation of A/8 between atoms in adjacent wells. 
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Appendix 

Using the Fourier series expansion 

w{x) = , (Al) 

n 

in Eq. H16|) we obtain the recursion relation 

A_in)W„.i + Ao(n)VF„ + A+(n)W„+i = 

S{p) [-B_^„,_i + Bo^n,0 + B+Sn,i] (A2) 
+ 2hka§ [B'_6n,-1 + B'oSnfi + B'Jr^,i] , 
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FIG. 8: Equilibrium energy in units of the recoil energy as a function of dimensionless intensity / for a = 10 and 5 = 0, 20, 35, 50. 



where 



l + ia{2n-l) 



1 + 2ina + (1 + ia{2n - 1))^ 



Ao{n) = 1 + 
A+in) 



a 



1 + 2ina 

,2 



l + ia{2n-l) 



1 + ia{2n+ 1) 



d2 + (1 + ia{2n - 1))2 d2 + (1 + ia{2n + 1))2_ 
1 + ia(2n + 1) 



1 + 2ina d2 + (1 + ia{2n + 1))2 



(A3a) 
(A3b) 
(A3c) 



1 d 

2 l-i2a d2 + (l-m)2 



a 



d 



2id2 + (l-za)2+d2 + (l + za)2/' (^'^^^ 
(T 1 d 

+ " 2 l + i2a d2 + (l + m)2 ' ^^^"^ 



" yi-z2ad2 + (l^- za)2 ' ^^^^^ 

^° = y{d^TTT^mF~d^TTTTmj^} ' ^"^^^^ 

, zcr 1 d 

+ " "yi + z2ad2 + (l + za)2 ' ^^^^^ 



15 



We are faced with solving the following equation: 
/ ■ 























■•. ^o(-2) A+{-2) 
A_(-l) Ao(-l) A+{-l) 








V 



A_(0) 












Ao(0) A+(0) 

^-(1) Ao(l) A+{1) 



where 



and j — ±,0. From Eq. 



A_(2) Ao{2) ■■. 

■■. 

Bj = S{p)Bj + 2hka—B'j , 

21, we see that for n < —1 and for n > 1 

A^in)Wn-i + Aoin)Wn + A+{n)Wn+i 



( 




\ 




( 




\ 









































So 






Wi 








B+ 



















\ 




) 




\ 




/ 



(A6) 



0. 



(A7) 



(A8) 



The final solution for the spatially averaged friction and diffusion coefficients depends only on Wq, W±i. However, to 
calculate these quantities all the other Ws must be evaluated. In practice, we truncate Eq. I|A6|I by setting W±n — 
and then compare the solution with that obtained by setting W±(^n+i) = 0; when these solutions differ by less than a 
fraction of a percent, we use the result to evaluate W±2/W±i, from which one can then calculate Wq, W±i. 
For n > 1, Eq. I|A6|I yields 



Wn 1 



A+in) Wn 



+1 



Ao(n) 



which can be written in the form 



A_(n)Mo(n) 



Setting n = 2 we obtain the continued fraction solution 

W2_ ^ 



1 



A+(«) 
Ao(n) W„ 



^-(2)Mo(2) 



A+(2) A_(3)Mo(3) 



1 



Similarly, for n < — 1 wc find 



W-2 



Ao{2) , A^(3) A_(4)Mo(4) 



A+(-2)Mo(-2) 



A-(-2) 
'4o(-2) i_ 



A+(-3)/Ao(-3) 

A_(-3) A^(-4)/Ao(-4T 
A0(-3) I---- 



(A9) 



(AlO) 



One can now use Eqs. 
one finds 



U, IjAlOjl to obtain equations for Wq, W±i in terms of C± and the BjS. Explicitly, 

(All) 









Wo 




Bo 









Ao(-l) + A_(-l)C_ A+i-l) 

^-(0) Ao(0) ^+(0) 

A_{1) Ao{l) + A+il)C+ 

The procedure is to obtain C+ and C_ according to the continued fraction solutions Eq. HA9|I and Eq. (|A10() and 
then find VF_i,o,i from Eq. IjAlip . 

Next we calculate 771,2, Ea. (|18a|) - (jl8b|) using Eqs. (|A1|) . IjAlip . (|14a|) as 



m 

m 



dS 

aoS + 2aihka— — h Woa2 + M^ifls + W^ia^ 
op 



boS + 2bihka 



dS_ 

dp 



WQb2 + Wib3 + W-ibi 



/2, 
/2, 



(Al2a) 
(A12b) 
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where 



ai 



dr e cos (dr) sin (ar) 



a(l-d2 + a2) 



/>oo 

/ dr e~'^ cos (dr) cos (ar) = — 
^0 (1 

POO 

02 = 2a / dr sin (dr) sin (ar) = — 
^0 (1 



(l + (d-a)2)(l + (d + a)2) ' 
(l + d2 + a2) 



+ (d-a)2)(l + (d + a)2) 
Aada 





POO 


as = 


id / dr e 




Jo 




POO 


04 = 


—ia I dr 




Jo 




ai , 


bi = 






POO 




2(7 / dr e 




Jo 




POO 


&3 = 


a dr 




Jo 




POO 


&4 = 


a dr 




Jo 



+ (d-a)2)(l + (d + a)2) 
iad 



d2 + (1 + ia)2 

—iad 
d2 + (1 - ia)2 



2crd(l +d2 -a2 



(l + (d-a)2)(l + (d + a)2) 
ad 



d2 + (l + ia)2 ' 
ad 



d2 + (l-m)2 ■ 

Since each Wi contains a term proportional to S and another term proportional to ^ , one has 

as 



W, = e,S + 2hkaQ 



dp ' 



and Eqs.(A12) can be written as 
= 



dS , 



m = 



S{ao + 6*002 + Oia^ + 9^ia^) + 2hka—{ai + (oa2 + Ci^s + C-i«4) 

S{bo + eob2 + 9ib3 + 0-164) + 2hka—ibi + Cob2 + C163 + C-164) 

op 



/2, 
/2. 



(Al3a) 
(A13b) 
(Al3c) 
(A13d) 
(Al3e) 

(Al4a) 
(A14b) 

(Al4c) 
(A14d) 
(Al4e) 

(A15) 

(Al6a) 
(A16b) 



Using Eqs. (AI6) with the Fokker-Plank equations H17() and 1)19(1 . one can identify the averaged force and diffusion 
coefficients appearing in Eqs. H2U|I with 



^/ = (flo + 000-2 + OlO-s + 0-104:), 

= {bo + dob2 + Oibs + O^ibi), 
Cst = (fli + Coa2 + CiOa + (-104)- 



(Al7a) 
(A17b) 
(Al7c) 
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